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ABSTRACT 



o 



Aims. We consider how tight binaries consisting of a super-massive black hole of mass M = 10^ - 10* Mq and a white 
dwarf in quasi-circular orbit can be formed in a globular cluster. We point out that a major fraction of white dwarfs 

■ tidally captured by the black hole may be destroyed by tidal inflation during ongoing tidal circularisation, and therefore 
$ZL(' the formation of tight binaries is inhibited. However some fraction may survive tidal circularisation through being spun 

up to high rotation rates. Then the rates of energy loss through gravitational wave emission induced by tidally excited 
pulsation modes and dissipation through non linear effects may compete with the rate of increase of pulsation energy 
due to dynamic tides. The semi-major axes of these white dwarfs are decreased by tidal interaction below a 'critical' 
^ ' value where dynamic tides decrease in effectiveness because pulsation modes retain phase coherence between successive 

' pericentre passages. 

Methods. We estimate the rate of formation of such circularising white dwarfs within a simple framework, modelling them 
^SJ ■ as n = 1.5 polytropes and assuming that results obtained from the tidal theory for slow rotators can be extrapolated 

^ ' to fast rotators. 

(^y) ' Results. We estimate the total capture rate as ~ ~ 2.5 • 10~^ Ml'^r^^'^yr'^ , where M4 = M/W^ Mq and ro.i is 

QQ I the radius of influence of the black hole expressed in units O.lpc. We find that the formation rate of tight pairs is 

, approximately 10 times smaller than the total capture rate, for typical parameters of the problem. This result is used 

■ to estimate the probability of detection of gravitational waves coming from such tight binaries by LISA. 

• ' Conclusions. We conclude that LISA may detect such binaries provided that the fraction of globular clusters containing 

black holes in the mass range of interest is substantial and that the dispersion velocity of the cluster stars near the 
, radius of influence of the black hole exceeds ~ 20km/ s. 

, Key words. Black hole physics- Globular Clusters- White dwarfs- Tidal Circularisation 
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1. Introduction play of two processes. These are the effect of distant two 
JH body gravitational encounters changing the orbital angular 
, There are some observational indications and theoret- momenta of the stars and an interaction associated with the 
ical suggestions that favour the presence of black holes in presence of the black hole which removes the orbital energy 
the mass range - 10 - 10 Mq in the centres of globular ^ star and is effective only when the orbital angular mo- 
clusters. The observational arguments supporting this hy- mentum is sufficiently small. This type of interaction may 
pothesis relate to kinematical phenomena observed in the result either through tidal interactions or by the emission of 
centres of some globular clusters (e.g. Gebhardt et al 2000, gravitational waves induced by the stellar orbital motion. 
Gebhardt, Rich, & Ho 2002) and the presence of X-ray 

sources not associated with the central nuclei in certain The efficiency of tidal interactions is determined by the 
galaxies (e.g. Fabbiano 1989, Matsumoto et al 2001, Ghosh ^atio of orbital pericentre distance to the tidal radius - the 
et al 2006 and references therein). There are also some the- letter being the distance from the black hole below which 
oretical models of the formation of such systems (e.g. Miller sigmficant disruption of the star through mass loss induced 
& Hamilton 2002). A review of the observational and the- by tides occurs. On the other hand, the efficiency of or- 
oretical aspects of this problem has been recently given by bital energy loss due to gravitational wave emission is de- 
van der Marel 2004 termined by the ratio of the orbital pericentre distance to 

In this Paper we assume that there is a black hole of gravitational radius of the black hole, 
mass M ^ 10'^ — 10'* Mq in a star cluster and estimate the Since the tidal radii corresponding to white dwarfs for 
rate of capture of white dwarfs of mass m by the black the range of black hole masses considered here are larger 
hole. The capture rate is always determined by the inter- than their gravitational radii, orbital energy is changed 
mainly through the action of dynamical tides. This is in 
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siding in galactic centres where emission of gravitational 
waves is more effective for changing the orbital energy of 
white dwarfs (e.g. Ivanov 2002, Preitag 2003 and references 
therein) . 

Since the relative contribution of two body gravitational 
encounters to the orbital evolution decreases very sharply 
for small angular momenta, a star with sufficiently small 
angular momentum loses orbital energy through tidal in- 
teraction while the orbital angular momentum remains ap- 
proximately unchanged. The latter occurs because the star 
cannot store a significant amount of angular momentum in 
comparison to the orbit. Accordingly, the orbital eccentric- 
ity decreases during this process which will be referred to 
as orbital circularisation. 

As a result of this process a tight quasi-circular orbit 
around a black hole may be formed. A white dwarf on such 
an orbit can emit gravitational radiation in the frequency 
band of order of lO^^Hz which is the most favoured for 
the planned LISA space borne gravitational wave antenna. 
In principal this is able to detect gravitational waves with 
dimensionless amplitude as small as 10~^^, for an observa- 
tional time of one year. That means that the presence of 
such a white dwarf can be detected from distances of order 
of ~ lO'^Mpc. Taking into account the fact that globular 
clusters form a very abundant population of cosmic objects 
such systems may contribute significantly to the budget of 
sources of gravitational radiation available for LISA. 

There is one principal obstacle inhibiting formation of 
a binary pair consisting of a black hole of mass M and a 
white dwarf of mass m in a tight orbit around it. Because 
it is produced through tidal interaction, its semi-major axis 
will be close to the tidal radius. There the ratio of orbital 
binding energy to gravitational energy of the white dwarf 
is of order of (M/m)^/^ » 1. Thus, for such an orbit to be 
produced, an amount of energy far exceeding the internal 
binding energy of the white dwarf must be removed from 
it. When tides are effective, orbital energy is transferred to 
pulsation modes and thence to the internal energy of the 
star. Therefore, without an effective energy loss mechanism, 
the white dwarf could be easily unbound and thus destroyed 
by the tidal input of energy or tidal inflation. 

Here we propose and discuss such a mechanism, which 
may, in principle, allow unbinding due to tidal energy input 
to be circumvented for a range of orbital parameters of 
the star. The operation of this mechanism depends on the 
interplay of several factors influencing the orbital evolution 
of star which we introduce below. 

The character of orbital evolution under the influence 
of tides is mainly determined by three factors: 1) the time 
scale for decay of the stellar pulsations excited by tidal 
interaction, 2) the orbital parameters of the star, and 3) 
the rotation of the star. 

Let us first consider the decay of pulsations in a non- 
rotating star. In this paper we assume that decay of stellar 
pulsations is caused both by the emission of gravitational 
waves that occurs because of the time-dependent density 
perturbations associated with them, and also by the dis- 
sipation of pulsation energy leading to its conversion into 
internal energy of the star. The latter decay channel is as- 
sumed to result through non-linear effects. Since its prop- 
erties are very poorly understood at the present time, we 
shall consider the corresponding dissipation time scale to 
be a free parameter. But we shall assume that it is larger 
than the orbital period of the star. However, this is not 



an essential assumption for orbital evolution at sufficiently 
large semi-major axes, see below. 

Initially the stellar orbit will be highly eccentric with 

a tidal interaction that excites stellar pulsations occurring 
impulsively at pericentre passage (e.g. Lai 1997, Ivanov & 
Papaloizou 2004 hereafter IP). For long decay time scales, 
pulsations will always be present in the star, and as a result 
of every periastron passage, a new perturbation excited by 
tides is added. 

When orbital semi-major axis is sufficiently large or the 
orbital period sufficiently long, it has been established that 
tidally induced changes to it cause the phase correlation 
between preexisting pulsations and freshly excited ones to 
be lost. Then, both the energy content of excited pulsa- 
tion modes and the orbital energy of the star evolve in 
a stochastic manner (eg. Kochanek 1992, Kosovichev & 
Novikov 1992, Mardling 1995). Under this evolution, the 
mode energy and the orbital binding energy of the star 
grow on average, with a part of the mode energy being 
transferred to the internal energy of the star. Therefore the 
stochastic exchange of energy between the orbit and stellar 
pulsations leads to a decrease of the orbital semi- major axis 
and period as well as an increase of the internal energy of 
the star. 

However, once the orbital period is sufficiently short, 
or the semi-major axis is below a critical value, pulsation 
modes can maintain phase coherence between successive 
pericentre passages, stochastic evolution ceases and dy- 
namic tides are expected to become less efficient. At this 
point, the tidal evolution rate is determined by the natural 
decay timescales of the pulsation modes. A steady pulsa- 
tion energy typical of that induced through one pericentre 
passage may be maintained, rather than the growth that 
occurs through the cumulative effect of mode energy in- 
puts proceeding over many pericentre passages when the 
semi-major axis is large. 

An important issue is whether the internal energy added 
to the star during the phase of stochastic evolution is 
enough to cause its destruction. For a non-rotating star 
the internal energy obtained by the star when the critical 
semi-major axis is reached is larger than the stellar gravita- 
tional energy. Therefore, such a star could be disrupted by 
tidal heating before the critical semi-major axis is reached. 

However, during orbital circularisation orbital angular 
momentum is also transferred to the star until an equilib- 
rium rotation rate is attained. This requires a rotation rate 
corresponding to corotation at periastron or faster. Thus 
the star can be spun up to high rotation rates. When mea- 
sured in terms of the amount of energy input per periastron 
passage, the efficiency of tidal interaction is minimised for 
such a rotating star. See for example figure 8 of Lai 1997 
and also IP. Accordingly the time scale of circularisation is 
increased. 

In this situation, when dissipation of the mode energy 
as a result of non linear effects is not effective, the time 
scale for transmitting orbital energy to the pulsation modes 
may be larger than the time scale for removal of the pulsa- 
tion energy through gravitational waves emitted because of 
the time-dependent perturbation of the star. In this situa- 
tion the white dwarf may reach the critical semi-major ax;is 
without internal dissipation causing tidal inflation because 
of cooling by emission of gravitational waves. In the oppo- 
site limit of effective mode dissipation and internal heat- 
ing, the rapidly rotating white dwarf may attain the crit- 
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Tidal capture into stoctiastic 
regime. Orbital evoiution of 
initiaiiy non rotating star is due 
to tides. Rate of evoiution 
decreases as star spins up. 



Outcomes of orbitai evoiution of a star 
spun up by tides 



Case a) 



Gravitationai waves iimit 
mode energy 



Case b) 



iVIode energy dissipates 
Star spins up further 



a1) 



Tides 
cause orbitai 
evoiution 



a2) 



Gravitationai 
waves cause 
orbitai evoiution 



b1) 



Tides cause 
orbitai evoiution 
Star is disrupted 



b2) 



Gravitationai 
waves cause 
orbitai evoiution 




Provided initial pericentre distance not too close the 
star leaves impulsive energy input (or stochastic 
regime) avoiding disruption 

> Orbit circularises 



Fig. 1. A schematic illustration of the possible orbital evo- 
lutionary paths of a white dwarf that is tidally captured 
by an intermediate mass black hole. The star is initially 
scattered into a highly elongated trajectory that evolves 
due to impulsive energy and angular momentum exchanges 
between the star and orbit that occur every pericentre pas- 
sage, the pericentre distance remaining very nearly fixed. 
As a result of these, the internal energy and angular mo- 
mentum increase. However, this process is slowed as the 
star spins up and under favourable conditions gravitational 
radiation can become important so that it controls the or- 
bital evolution and/or damps the excited pulsation modes. 
The latter case a) leading to cases al) and a2) occurs where 
other non linear processes are ineffective at dissipating the 
pulsations. Case b) leading to cases bl) and b2) occurs 
when such processes are more effective. For both cases a) 
and b) the star may survive the orbital evolution if the tidal 
capture occurs at sufficiently large pericentre distance. 



ical semi-major axis through the emission of gravitational 
waves induced by orbital motion rather than through tidal 
interaction. 

Thus taking into account the reduction in effectiveness 
of the tidal energy transfer brought about by the effect of 
stellar rotation, we find that for sufficiently large orbital an- 
gular momenta the semi-major axis can be decreased below 
the critical one without significant heating of the star. 

After the critical semi-major axis is reached, again, be- 
cause of the reduced efficiency of the tidal interaction, or- 
bital evolution is governed by the emission of gravitational 
waves determined by the orbital motion of the white dwarf. 
This process can further reduce the orbital semi-major axis 
and lead to formation of a tight quasi-circular orbit. 



We formulate the criterion for white dwarf 'survival' 
treating dynamic tides within the framework of the sim- 
plest possible model of tidal interactions. We assume that 
the internal structure of a white dwarf is the same as that 
of a n = 1.5 polytrope. We also assume that the results 
obtained from the theory of dynamic tides in slowly rotat- 
ing stars can be extrapolated to high rotation rates for the 
purpose of making approximate estimates. Based on these 
assumptions we estimate the formation rate of tight pairs 
which turns out to be an order of magnitude smaller than 
the total rate of tidal capture, for typical parameters of 
the problem. This allows us to obtain an estimate of prob- 
ability of detection of such sources of gravitational waves 
by LISA. We conclude that LISA could, in principal, de- 
tect such a source provided that there is a significant frac- 
tion of globular clusters containing black holes with masses 
~ 10^ — IO^Mq, and with stellar velocity dispersions in 
their innermost regions exceeding ~ 20km/ s. 

There are a number of different possibilities and branch- 
ings associated with the orbital evolution subsequent to 
tidal capture. Each of these requires consideration of sev- 
eral physical processes and is described by a number of 
algebraic expressions. In order to clarify the situation, we 
provide a more transparent summary of the proposed paths 
to a circularised orbit with disruption of the star avoided. 
In Figure[T]we give a diagrammatic illustration of the possi- 
ble orbital evolutionary paths that can be taken by a white 
dwarf subsequent to tidal capture by an intermediate mass 
black hole. Initially, impulsive energy and angular momen- 
tum exchanges between the orbit and star that occur every 
pericentre passage spin it up and excite modes of oscilla- 
tion. As a result the rate of tidal evolution of the orbit de- 
creases such that gravitational radiation can become more 
important. This may also be important for damping the 
oscillation modes. If the tidal capture occurs at sufficiently 
large pericentre distance, the importance of gravitational 
radiation during the orbital and pulsation mode evolution 
may allow the star to survive the regime of impulsive en- 
ergy input at pericentre passage and become circularised 
without disruption and be a potential LISA source. 

The plan of the paper, in which the above phenomena 
are considered in more detail, is as follows. In Section [5] 
we describe the model of the stellar cluster we use and the 
effect of distant two body encounters. In Section [3] we go 
on to discuss the tidal interaction of a white dwarf with the 
central black hole and pulsation mode excitation through 
dynamic tides. Then in Section |4] we consider the way in 
which gravitational wave emission/non linear effects may 
lead to pulsation mode amplitude limitation. In Section [5] 
we go on to explore the conditions for safe circularisation or 
tidal disruption of the white dwarf. In Section [S] we discuss 
the rate of tidal capture of white dwarfs by the black hole 
and estimate the fraction that can circularise safely, with 
technical details being relegated to an Appendix. In Section 
[7] we estimate the probability of detection of this source of 
gravitational radiation. Finally in Section [8] we discuss our 
results. 



2. The stellar system 
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2.1. Density distribution 

In order to discuss the rate of orbital circularisation for 
white dwarfs it is first necessary to consider the properties 
of the star cluster in which they are situated. This is be- 
cause dynamical interaction with the stars in the system 
is responsible for the production of orbits for which tidal 
interaction with the central black hole becomes significant. 

The radius of influence of the central black hole, r^, is 
defined as the radius within which the total mass of the 
stars is the same as the mass of the central black hole: 
Mst{ra) = M. As has been noted by many authors, when 
r < Ta, a 'cusp' in the stellar distribution of stars is formed 
with a density distribution such that 



p cx r 



-(3/2+p) 



(1) 



where the parameter p is typically small. Two values that 
have been suggested through theoretical considerations are 
p = (Young, 1980) and p = 1/4 (Gurevich, 1964, BahcaU 
& Wolf 1976). Recently simulations by Baumgardt, Makino 
& Ebisuzaki 2004a, and Baumgardt, Makino & Ebisuzaki 
2004b have indicated that for the case of globular clusters 
having a stellar population with a single mass, p ~ 1/4, 
while for systems with a realistic multi-mass stellar popu- 
lation, p w 0.05, with p « 1/4 for the most massive stars. 

In this paper we make only approximate estimates of 
the tidal circularisation rate. Thus for simplicity, we have 
adopted p = for all components of the stellar system. 
We note that if a steeper cusp profile were to be adopted, 
the number of tidally captured stars at any particular time 
would increase by a numerical factor of order unity. 

2.2. Pliase space distribution 

Making the usual assumption that the phase space distri- 
bution function for 'typical' stars depends on the binding 
energy per unit mass only and the distribution of white 
dwarfs is the same as that of typical stars, it can be easily 
shown from ([T]) with p = that the distribution function 
for the white dwarfs over binding energy per unit mass and 
specific angular momentum, N^j^{E, J), is 



\mriJ 



(2) 



where M and m are the masses of the black hole and of 
a typical star, and we assume that m « Mq. 6 0.1 is 
the number fraction of white dwarfs in the cusp, E and 
J are, respectively, the binding energy per unit mass and 
the specific angular momentum of a star. We adopt the 
convention that the binding energy has positive values for 
the gravitationally bound stars considered in this paper. 



PorbiE) 



TT GM 



(3) 



is the orbital period of a star with binding energy per unit 
mass E. 



2.3. Loss cone 

The assumption of isotropy of the phase space distribu- 
tion function is applicable only for estimates of quantities 
determined by the bulk of the stars, such as e.g. quan- 
tities characterising distant gravitational interactions, see 



equation ([5]) below. It becomes invalid for stars with suf- 
ficiently low orbital angular momentum which can either 
be tidally disrupted or directly captured by the black hole. 
The presence of the black hole thus leads to the formation 
of a 'loss cone' such that equation ^ over-estimates the 
number of stars having low angular momenta. Assuming 
that stars with sufficiently small specific angular momenta 
J < Jcone are absent, the presence of the loss cone may be 
easily accounted for by a correction factor (e.g. Lightman 
& Shapiro 1977) such that the distribution function of the 
white dwarfs is modified to become 



wd 



(4) 



where 9 = J/ Jcone, Ai = \n{Jctrc{E)/ Jcone), and Jcirc is 
the specific angular momentum corresponding to a circular 
orbit with a given value of the binding energy per unit mass 
E. Typically, Ai - 8. 

2.4. Distant encounters 

A star orbiting around the black hole changes its energy 
and angular momentum due to distant gravitational in- 
teractions with other stars in the cusp. These result in a 
random walk in phase space. The stochastic change of the 
specific angular momentum per orbital period may be char- 
acterised by its dispersion, (Aj2)^, which can be written in 
the form (e.g., Ivanov, 2002): 



-5/2 



(5) 



where q = m/M is the mass ratio, A2 « ln(0.5g ^) and 
we use the dimensionless energy variable £ defined through 
£ = Era/{GM)\^ 



3. Tidal interaction of a white dwarf with the 
central black hole 



3.1. White dwarf model 

In order to make order of magnitude estimates we use the 
simplest possible model for the white dwarf. We assume 
zero temperature so that the pressure is due to completely 
degenerate electrons. In this case the equation of state is 
baratropic and gravity or g modes are absent. Furthermore 
we approximate the structure of the white dwarf as that of 
an — 1.5 polytrope. In addition to the mass m and radius 
Rwd we introduce the two parameters /x = m/MQ and A = 
^u,d/(0.01i?o) i?„,d/(7-10*cm). The effectiveness of tidal 
interactions depends significantly on the mean density of 
the white dwarf. To characterise low and high density cases 
we consider two sets of values for /i and A. The 'low density' 
case will be characterised by parameters appropriate for a 



^ In what follows it is convenient to represent different en- 
ergy scales in dimensionless form. There are two natural units 
of energy per unit mass in our problem: the first associated the 
stellar cusp - GM/va and the second associated with the in- 
ternal structure of star itself - Gm/r^d. For the dimensionless 
quantities expressed in terms of the first unit we will use the 
calligraphic style while the quantities expressed in terms of the 
second unit will be denoted by tilde. 
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'typical' white dwarf for which = 0.6 and A = 1.4, and 
for the 'high density' case we use parameters corresponding 
to Sirius B for which w 1 and A « 0.84. 



3.2. Tidal perturbations and orbital circularisation 

For a baratropic star only fundamental (/), pressure (p) and 
inertial modes can be excited by tidal interaction. Since p- 
modcs have large cigcn-frequencies and inertial modes are 
significant only for rather large orbital angular momenta 
(Papaloizou & Ivanov 2005), their respective contributions 
to the tidal exchange of energy and angular momentum are 
small compared to the contribution of the /-modes for the 
range of the orbital parameters of interest, so they are not 
considered further. 

To make our estimates, we use the linear theory of tidal 
perturbations induced in a slowly rotating star and consider 
only the leading order quadrupole response. The assump- 
tion of slow rotation is certainly valid in the early stages of 
orbital evolution due to tides. In this regime orbital changes 
induced by distant gravitational interactions and tidal ef- 
fects compete with each other. Then we may assume that 
the star is non-rotating when calculating the rate of tidal 
circularisation. 

However, the star can be spun up because of the effect of 
tides during orbital circularisation, achieving a considerable 
rotation rate during the later stages. As we will see below 
this may have important implications for the outcome of 
the orbital evolution of the star. In particular, it may play 
a role in determining whether is it possible or not for the 
star to survive in a tight quasi-circular orbit. Therefore, it 
is directly related to the problem of formation of sources 
of gravitational radiation. Unfortunately, a theory of tidal 
perturbations of a fast rotating star is practically unde- 
veloped in the present time. Accordingly, we shall assume 
that the results obtained from the theory of slowly rotating 
stars can be extrapolated to the case of fast rotation for 
the purposes of making order of magnitude estimates. 

During the early stages of tidal circularisation, the orbit 
of the star is highly eccentric. In this case tidal interactions 
are only important near pericentre. Because of this, in order 
to calculate the energy and angular momentum exchange 
between the orbit and the modes of pulsation of the star 
it is possible to consider it as undergoing a sequence of 
successive flybys of the black hole each of which produces an 
impulsive change. This formulation of the problem allows 
us to use the theory of tidal excitation and dissipation of 
the fundamental modes in a baratropic star moving on a 
highly eccentric or parabolic orbit developed elsewhere (see 
e.g. Press & Teukolsky 1977 hereafter PT, Lai 1997, IP). 
This theory is based on consideration of a single pericentre 
passage where an amount of energy AE and an amount 
of angular momentum AL is transferred from the orbit to 
the star. These quantities depend on the amplitudes and 
phases of the pulsation modes before the passage and on 
the radius of pericentre, rp. 

Let us temporarily assume that the amplitude of the 
pulsation mode before pericentre passage is negligible and 
the stellar model and the angular velocity of the star are 
fixed. In this case, it can be shown that the quantities AE 



and AL are determined by the value of dimensionless pa- 
rameter (PT) 

where rp is the pericentre distance. These quantities can 
be naturally considered as being made up of two parts: 1) 
a quasi-static part determined by irreversible loss of total 
energy of the system during the pericentre passage and 2) 
a dynamic part (due to the so-called dynamic tide) deter- 
mined by excitation of the normal modes of the star. 

3.3. Quasi-static tides 

The quasi-static part of the energy transfer AEq_g, can 
be written in the form: AEq-g « 'yE^./rf', where = 
Gm^ I R^d is some 'reference' energy associated with the 
star and the dimensionless coefficient 7 is determined by 
the dissipative processes operating in the star. An explicit 
expression for 7 can be found in IP. During the flyby energy 
loss occurs as a result of viscous dissipation in the bulk of 
the white dwarf and in its convcctive envelope and also due 
to emission of gravitational waves. Since the coefHcient of 
dynamic viscosity in the bulk of the white dwarf appears 
to be very small (^ 10^ — 10^ in cgs units, e.g. Chugimov & 
Yakovlev 2005 and references therein) and convcctive en- 
velopes of white dwarfs also have a rather small estimated 
turbulent viscosity and a very small relative mass, their re- 
spective contributions to the coefficient 7 are quite small. 
The main contribution appears to come from the emission 
of gravitational waves generated by time dependent per- 
turbation of the white dwarf (see Osaki & Hansen 1973, 
hereafter OH). However, a simple estimate shows that the 
value of 7 determined by this process leads to a value of 
AEq^s much smaller than that associated with dynamic 
tides, for relevant values of 77. Therefore, the overall con- 
tribution of quasi-static tides appears to be negligible and 
will not be considered further. 

3.4. Dynamic tides 

Now let us consider dynamic tides associated with the fun- 
damental quadrupole mode of pulsation. As was shown by 
IP for the case of a slowly rotating star, the quadrupole 
mode propagating in the direction of orbital motion (the so- 
called prograde mode) determines the exchange of energy 
between the orbit and the star provided that the angular 
frequency of the star, O^., is smaller that its 'equilibrium' 
value Vieq « (2 -I- ln?7)f7,/77, where 1]» = Gm/Rl^^ is a 
'reference' frequency associated with the star. As it will be 
shown later, in our problem a typical value of ~ 4, and 
the equilibrium angular frequency fieg ~ O.Sfi* is formally 
larger than the angular frequency at rotational break-up of 
the star, 17^,. w 0.551,. Although the theory leading to this 
conclusion is not valid at such high rotation rates, it seems 
reasonable to assume that the energy exchange is mainly 
determined by the prograde mode unless a rate of rotation 
very close to Q.hr is reached. Considering excitation of only 
this mode and assuming that the star is not perturbed be- 
fore pericentre passage, it easy to see from results given 
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in IP that the energy gained by the star per unit of mass, 
AEt can be expressed as 



R 



■wd 



(7) 



where we use the fact that the eigenfrequency of the 
quadrupole / mode for a non-rotating n=1.5 polytrope 
w/ « 1.445r2, and the value of dimensionless 'overlap' inte- 
gral Q/ ~ 0.5 calculated in, e.g., Lee & Ostriker 1986. The 
quantity 0(?7) characterises the strength of the tidal inter- 
action for a non-rotating star and ^' determines correction 
to the energy exchange due to rotation, being unity for a 
non rotating star (see below). Explicitly, we have 



0(?7) =?7'iexp(2.74(r/-l)), 



(8) 



where rj is given by equation When ^ 1 — 2, the white 
dwarf is either tidally disrupted or stripped, and the linear 
theory of dynamic tides cannot be used (see e.g. Frolov et 
al. 1994; Ivanov & Novikov 2001 and references therein). 
As it will be shown later, typically, 77 ~ 4 — 5 for white 
dwarfs undergoing tidal circularisation. Since the energy 
gain decreases exponentially with rj the use of the linear 
theory seems to be justified in our case. 

3.5. Direct capture or tidal disruption 

The characteristic pericentre distance below which tidal 
disruption is expected is determined as the distance where 
the tidal force acting on an element of the star near its sur- 
face is comparable to the gravitational force due to the star 
itself. Thus it is given by 



M 



1/3 



= ( — ) Rwd ~ 1.5 • 10^" Xfi-^^^A'll^^ cm 



(9) 



where M4 = A//(10^A^q). Using this we can rewrite equa- 
tion ^ in another form 



J 



(10) 



where Jt = \/2GMrT is the specific orbital angular mo- 
mentum of an orbit with pericentre distance equal to rx- 
From equation (jlOp it follows that when J w const rj oa 
is also approximately conserved. 

The black hole mainly tidally disrupts the white dwarfs 
provided that its mass is sufficiently small. To see this we 
note that the specific orbital angular momentum of stars 
directly captured by the black hole must be smaller than 
Jcap = 4:GM/c. Define the 'critical' mass Merit by the con- 
dition that Jcap = Jt when M = Merit- One obtains 



Mr. 



1.4 • lO-^A^/ V^^^M 



(11) 



When M < Merit, the specific angular momentum Jt ex- 
ceeds Jcap so it determines the boundary of the loss cone. 
Thus Jcone ~ Jt and, accordingly, 77 « 0^ ( recall that 6 
is defined as 6 = J/ Jcone)- 



3.6. Effect of stellar rotation 

As indicated above, the rotation rate of the star has an 
important effect in determining the strength of the tidal 



interaction. From the results given in IP it follows that the 
corresponding correction factor ^E* entering in equation ([7]) 
has the form 



4V2~ 

exp "r??, 



where fi^ = ^r/^*, ^r 
and we recall that = 



(12) 

is the angular velocity of the star 
characterises a typical fre- 



quency of stellar pulsation associated with the fundamental 
mode. It is important to note that equation p2p has been 
formally derived assuming that VLr <^ 1. However, we are 
going to use it later on to make our estimates implying that 
it is approximately valid when f2j. < f2t,r ~ 0.5$!*. Note 
that this assumption may be justified in part by results 
of calculations of Managan (1986) who considered normal 
modes of a fast rotating n — l.b polytrope and of Clement 
1989 who considered those of a fast rotating ISMq star. 
In both cases it was found that the eigen frequency of the 
prograde / mode has an approximate linear dependence on 
the frequency of rotation until an angular frequency very 
close to ^br is reached. The linear dependence of the eigen 
frequency on the frequency of rotation appears to be one 
of the most important factors leading to equation ([T^ . 

3. 7. The effect of many pericentre passages and the 
criterion for stochastic instability 

Let us consider a sequence of many pericentre flybys as- 
suming that the energy per unit mass contained in modes 
of oscillation, Em, does not decay significantly between suc- 
cessive pericentre passages. In this case there is interference 
between the preexisting wave perturbation in the star and 
the wave excited in the vicinity of pericentre as a result 
of tidal perturbation. In this case, simple addition of ([7]) 
to Em after the pericentre passage to obtain the new Em 
needs further justification. In this context it is important 
that the orbital period Porb is changing with time as a result 
of tidal interaction and loss of orbital energy due to emis- 
sion of gravitational waves 0, and therefore the the phase 
change of the oscillating mode generated between two suc- 
cessive pericentre passages, $ = LOfPorb, is also changing 
with time. When the change of $ during one orbital pe- 
riod, 5$ — ujfSPorb is larger than 27r, the mode phase at 
the time of the next pericentre passage is essentially a ran- 
dom quantity, and accordingly, the mode energy per unit 
mass Em evolves as a stochastic variable. In this case we 
can add expression ([7]) to Em at each pericentre passage on 
average, assuming that all resulting expressions are valid in 
some statistical sense being averaged over many realisations 
of the stochastic process. The process of stochastic evolu- 
tion of the mode energy will be refereed to as 'stochastic 
instability'. 



^ Note that it is important to distinguish between the effect 
of emission of gravitational waves by perturbations of the white 
dwarf and the effect of emission of gravitational waves due to the 
orbital motion of the star. The former effect is important as a 
major process of mode energy loss while the latter significantly 
influences the orbital evolution and may provide a source of 
gravitational waves during the late stages of orbital evolution, 
see below. 
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A condition for the stochastic instability to operate can 
be readily found from the condition lo f\dPorb / dE\/^E > 2tt 
and the expression ([3]) for the orbital period in the form 

a>a.t^ J^^^ - A^-^/^Af^ (A^)-^/5l0"cm, (13) 
where a is the orbit semi-major axis, 

Rwd 



oscillating mode and the resulting decay time can be ex- 
pressed as 



tde 



1.5 • 10^X^n--'yr. 



(15) 



AE = AEt + AE, 



GW 



Gm 



{AEt + AE, 



GW) 



(14) 



is the dimensionless change of the orbital energy per unit 
mass, AE, in one orbital period, and AEqw is the change 
of the orbital energy per unit of mass during one orbital 
period due to emission of gravitational waves, see equation 
l[28|) below for an explicit expression. As follows from equa- 
tion (HH) the condition for the stochastic instability to set 
in involves two contributions to the change of the orbital 
energy: AEt associated with tides and AEqw determined 
by gravitational waves emitted as a result of orbital mo- 
tion. This differs from the case of tidally interacting stellar 
binaries where only the tidal contribution is normally con- 
sidered. In the case of stellar binaries, when the stochastic 
instability operates, both the orbital binding energy and 
the energy of the mode of pulsation evolve in a stochastic 
manner (e.g. Mardling 1995). The fact that AEqw enters 
into equation leads to the possibility of the stochastic 
instability operating in a situation where orbital evolution 
is mainly determined by emission of gravitational waves, i.e. 
AEgw ^ AEt, see Section 5 and the Appendix. In this 
limit, our numerical investigations indicate that the orbital 
binding energy can increase in a way that appears to be 
fairly regular while the energy of the pulsation mode still 
tends to increase in an erratic manner. On the other hand, 
the mode energy ceases to increase when a <^ ast and 
AEt / AEqw >^ 0.05. We have checked that the condition 
AEt / AEgw >~ 0.05 is fulfilled for parameters typical for 
our problem. 

It is interesting to note that the condition (|13p can be 
also derived from quite a different approach based on the 
treatment of the tidal interaction as occurring through res- 
onances between the mode frequency ujf and the chang- 
ing orbital frequency fiorb of different order n, such that 
nriorb = ^f- The condition of overlap of these resonances 
given that florb is changing is Sn = uj f\dil,~^fj/ dE\AE > 1. 
This is equivalent to the condition (fTF 



4. Gravitational wave emission/non-linear effects, 
characteristic time scales and mode amplitude 
limitation 

As we have mentioned above the most important linear de- 
cay channel of the fundamental mode appears to be through 
emission of gravitational waves. The corresponding decay 
time of the fundamental mode has been calculated by OH 
for two models of iron white dwarfs. The results can be 
extended to the range of masses and radii appropriate for 
'standard' CO white dwarfs with help of the expression for 
the gravitational wave luminosity, Lqwi produced by the 



^ A more accurate expression has been derived in IP. See also 
this paper for the references on numerical and analytical studies 
of the stochastic instability. 



Note we define the decay time as tdec = Em/Law which 
gives a value of t^ec which is a factor of two smaller com- 
pared to that of OH. 

In this situation the mode energy can attain a large 
value as a result of the cumulative effect of multiple tidal 
interactions occurring at pericentre passage. At high ampli- 
tudes non linear effects may lead to additional dissipation 
and mode decay. We characterise such effects by adopting a 
mode decay time scale tni which is a function of the mode 
energy. The non linear theory for these pulsation modes 
is not yet adequately developed to provide a form for t„;. 
Therefore, we regard this as a free parameter and consider 
two limiting cases: a) The time tni is assumed to be large 
compared to idee, being taken to be comparable to te«, the 
actual time scale of orbital evolution which is the smaller of 
the times required to change the orbit significantly through 
tidal interaction or to make it evolve as a result of the emis- 
sion of gravitational waves induced by that orbital motion, 
see below equation (P^ . Assuming that all quantities char- 
acterising the mode of pulsation and the star evolve on 
the time scale of order of tg^, we estimate below that the 
condition tni ~ tev is equivalent to the condition that the 
rotational angular momentum of the star is the same order 
of magnitude as the angular momentum corresponding to 
the mode, see equation (|18p below. Thus, in this situation 
the degree of non-linearity adjusts in a way that does not 
allow the mode angular momentum to significantly exceed 
that in the stellar rotation, b) The time tni is assumed to 
be short compared to both tdec and the orbital evolution 
time scale tev However, in both cases we shall assume that 
tni is much larger than the orbital period of the tidally in- 
teracting star. But note that this is not strictly necessary 
during the first phase of evolution when the semi-major axis 
is large ( see section 14.11 below) 



4.1. Balance of mode energy and angular momentum input 
due to tides and losses due to gravitational 
radiation/non linear effects 

Here we estimate the amount of mode energy that is stored 
in the star when there is a balance between the build-up of 
mode energy due to the stochastic instability and decay due 
to the emission of gravitational waves or non linear effects. 
Provided this balance can occur with a stellar mode angular 
momentum content that is small enough to avoid break up 
of the star, it is possible for it to avoid disruption by having 
to absorb a large amount of released orbital binding energy. 
In a stationary state the mode energy per unit of mass can 
be estimated from the balance equation 



E„ 



1 



1 

tni 



AEj 



(16) 



The mode specific angular momentum, L^, is related to 
Eni by the standard equation: 



,Em 



(17) 



where we take into account that tides mainly excite the 
quadrupole mode with the azimuthal mode number m — 2. 
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The rate of transfer of angular momentum to the star 
is related to the rate of dissipation of mode energy (eg. 
Goldreich & Nicholson 1989) and thus determined by tni- 
The rate of transfer of specific angular momentum is L* = 
Lm/tni and accordingly the total specific angular momen- 
tum transferred to the star in a time interval tev is esti- 
mated to be 



tni 



(18) 



Obviously, should not exceed the specific angular mo- 
mentum of the star in a state of uniform rotation with an- 
gular frequency corresponding to rotational break-up. 

Estimating the moment of inertia of the star per unit 
of mass as / w 0-2Rla and using the standard relation 
L* = IQ^,Qr, from equations (fT6|) . (fT7|l and (fT8|l we obtain 
the associated dimensionless angular velocity of the star to 
be given by 



7- 



tde 



(19) 



{Ud + tdec) V *T 

where we introduce the characteristic time for tidal build- 
up of the mode energy 

E 



tT = 



AEt 



-P, 



orb • 



(20) 



Note that tx so defined can formally be smaller than Porb- 
It is obvious that in that case it is the orbital period 
which plays the role of the characteristic tidal time scale. 
Substituting equation ([7]) in (^0)1 we can express the tidal 
time scale in the form 



tT 



(21) 



where ri = ra/{lpc). Note that since {(j)'^) ^ 1 typically 
we have tx ^ ^yr. 

In general the time scale of orbital evolution is given by 

tei, « min(iT, ^Gw), (22) 

where taw is the characteristic time of orbital evolution 
due to gravitational wave emission explicitly given below, 
see equation (|29p . As indicated below, for a white dwarf 
spun up through tidal interaction, gravitational waves emit- 
ted as a consequence of orbital motion are more effec- 
tive than tides at inducing orbital evolution. In this case 

tev ~ taw ^ tr- 

Recalling that the break up angular velocity Qi,r ~ 
0.5f2*, and assuming that the criterion fir < Qbr provides 
the condition for the white dwarf to survive the circularisa- 
tion process, we obtain a limitation on the orbital binding 
energy per unit mass of the white dwarf, during the phase 
of stochastic instability, to be given by 



E < Ed^s = 7-10- 



> {tni + td ec) 

^dec 



tr 



Gm 

Ryjd 



(23) 



This equation leads to a constraint that the orbital binding 
energy of the white dwarf should not be too high while 
mode energy input as implied by the existence of stochastic 
instability is occurring. The condition (|23p can be simplified 
when tni is either large or small. For case a) when tni ~ t^v 
the condition can be written as 



dis 



= 7-10 



_2 m rg tr 
M R^d td ec 



(24) 



where we recall that the dimensionless energy £ is re- 
lated to the orbital binding energy per unit mass E as 
£ = {'ra/GM)E. This condition can be reformulated as 
setting an upper bound on the allowed mode energy per 
unit mass as 



E,n < 7 • 10" 



, Gm 

Rwd 



(25) 



Note that in case a) the angular momentum transferred to 
stellar rotation is approximately equal to that associated 
with the mode of pulsation. Therefore the condition tni ~ 
tgy can be regarded as a condition for equipartition of the 
angular momentum associated with the pulsation mode and 
that associated with the stellar rotation. 

In the opposite limit corresponding to case b) we have 
tni ^ tdec, tev Then equation gives 



£ < £d^s = 7-10- 



tr 



M Ry,d t e 



(26) 



In this case the energy put into internal energy of the star, 
Eint , through dissipation is much larger than the mode en- 
ergy, Errn and we can express (pS)) as leading to an upper 
bound for Eint in a form analogous to that given by (|25p 



E,nt < 7 • 10- 



Gr 



R 



■wd 



(27) 



We remark that this condition, obtained from consid- 
erations of the stellar rotation, is equivalent to specifying 
a limit on the allowed increase of the internal energy. As 
almost all of the orbital binding energy is converted to inter- 
nal energy of the star in case b) , the same constraint would 
occur if the initial evolutionary phase involving stochastic 
instability were replaced by one in which the energy input 
into oscillations during each pericentre passage was dissi- 
pated before the next one. Thus, as indicated in section |3] 
above, the condition that tni should be less than the orbital 
period is not strictly necessary during the initial evolution- 
ary phase following tidal capture (see also IP). 

Note too that we have neglected the change of the white 
dwarf radius due to tidal heating. Assuming this change 
is modest, we can take it into account by considering the 
white dwarf to have a slightly smaller mean density when 
discussing conditions for safe circularisation, see the next 
Section. 



4.2. Effect of gravitational radiation on the orbit 

The change of orbital energy per unit mass due to emis- 
sion of gravitational waves generated by orbital motion of 
the star can be obtained from results given by Peters 1964. 
In the case of a highly eccentric orbit, the change per or- 
bital period, which is mainly induced at pericentre may be 
written in the form 



AEaw = 3.8 • 10" 



J 



GM 



io-V'/'a-5/2//6m4/3, 



(28) 



and the corresponding evolution time can be expressed in 
terms of as 



taw 



AErr 



AE, 



tT. 



(29) 
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Since the corresponding quantity induced by tides, AEt, 
decreases with i] exponentially and AEqw as an inverse 
power of 77, there is a certain value of 77 where tr = taw- At 
larger values of 77 the orbital evolution is mainly determined 
by emission of gravitational waves. For non rotating stars 
(i.e. setting ^ = 1) this value of ry is approximately equal to 
~ 5.3 for a 'high density' white dwarf and ~ 5.6 for a 'low 
density' white dwarf. As we have mentioned above and will 
discuss below, a typical value of rj for the stars undergoing 
tidal circularisation is of the order of ~ 4 and always smaller 
than 5. Taking into account that tx increases very rapidly 
with this parameter, we can conclude that tidal interactions 
play the most significant role in the formation of a flow of 
stars undergoing orbital circularisation. 

On the other hand, the star is spun up by tides dur- 
ing ongoing orbital circularisation and It is increased as 
a result. Therefore, at a later stage, the orbital evolution 
may be governed by emission of gravitational waves when 
the orbital binding energy is sufficiently large, depending 
on parameters of the white dwarf and of the star cluster. 
Moreover, when the binding energy exceeds the value cor- 
responding to the onset of stochastic instability, the tidal 
response becomes quasi periodic and ineffective so that the 
orbital evolution is solely determined by the emission of 
gravitational waves. 

5. Conditions for safe circularisation or tidal 
disruption 

Before going on to estimate of the rate of tidal capture and 
circularisation we use the results of the previous sections 
to examine more closely the conditions required for mode 
amplitude limitation due to gravitational wave emission to 
enable safe circularisation with the possibility of tidal heat- 
ing, inflation and disruption avoided. 

During the early evolutionary phases, subsequent to 
tidal capture, the star's orbit is highly eccentric, and the 
orbital angular momentum is approximately conserved dur- 
ing the orbital evolution while the semi-major axis/binding 
energy as well as the angular velocity of the star change 
with time (see eg. IP). When tides dominate over the emis- 
sion of gravitational waves induced by orbital motion, the 
evolution time scale of the semi-major axis is given by equa- 
tion (pn|) and we have t^y « tT- When gravitational radia- 
tion controls the orbital evolution, which happens at large 
enough orbital binding energy, we have t^v ~ icw where 
taw is given by equation P5)) . For a given value of the or- 
bital angular momentum, the orbital evolution time t^-y de- 
pends on the orbital energy E, and thus, on the star's semi- 
major axis a. This dependence can be found from equations 
pg)) . pi)) . dUl), (123 and (1211). For illustrative purposes we 
plot this dependence as well as that of other characteris- 
tic time-scales in Figures 1 and 2 for a 'low density' white 
dwarf orbiting around a 10** M0 black hole. The cusp size is 
taken to be Ipc. All dependencies are calculated in assump- 
tion that the stochastic instability operates for the shown 
range of binding energies. 

In Figure 1 the characteristic time-scales corresponding 
to a star with 77 = 5 are shown. The horizontal line indi- 
cates the decay time scale of the stellar pulsations due to 
the emission of gravitational waves, tdec^ given by equation 
p3)) . The inclined dotted line indicates the dependence of 
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(r,/GM)E 

Fig. 2. Characteristic time scales years as functions of the 
dimensionless orbital energy £. Plotted here are tdec, the 
decay time of the pulsation mode due to emission of grav- 
itational waves associated with the stellar pulsation, te„, 
the time scale for evolution of the orbit. It-, the time scale 
of build-up of the mode energy and tow , the time scale 
of orbital evolution due to emission of gravitational waves 
induced by orbital motion. The horizontal dotted line 1 rep- 
resents idee- The two solid curves 2a and 2b represent t^y 
calculated assuming that the stellar rotation frequency is 
given by equation (jl9p . The curve 2a taking on larger val- 
ues at small values of E represents the case tjix ~ ^ei) , and 
the other curve represents the case tni <C t^v The inclined 
dotted line 3 represents which becomes equal to t^v 
at the larger binding energies where gravitational radiation 
controls the orbital evolution. The dashed and dot dashed 
lines 4 and 5 represent It calculated for $1^ = f2f,r and 
= 0, respectively. The dot dot dashed curves 6a and 6b 
represent It calculated for the star rotating with angular 
frequency given by equation for values of binding ener- 
gies larger than that corresponding to the intersection of t^v 
and taw- For smaller values of energies we have It = ie-u- 
Note the slower evolution due to tides for the faster rota- 
tion rate. This results in gravitational radiation ultimately 
controlling the orbital evolution. All curves are calculated 
for II = 0.6, A = 1.4, 77 = 5, ri = 1 and M4 = 1. 
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Fig. 3. Same as Figure 1 but for t] — 4. Note that in this 
case we have tev = tx for the whole range of binding ener- 
gies shown. 
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the characteristic time taw on £ and the inclined dashed 
and dot-dashed lines show the dependence of tx on £ for 
fir = ^Ibr = 0.517* (dashed line) and fir — (dot-dashed 
line). We recall that the timescale tr applies to the evo- 
lutionary phase just after tidal capture which can exhibit 
stochastic instability. It may not be applicable to the later 
stages, when £ is large, and the conditions for additive im- 
pulsive energy inputs at pericentre passage are not satisfied, 
see discussion below. 

For a given value of the dot-dashed line corresponding 
to a non rotating star, gives the smallest possible value 
of It = tmin- For fixed £, stars with values It exceeding 
that given by the dashed line, would have angular velocities 
exceeding Qbr, and so would be disrupted. Therefore the 
dashed curve gives the largest allowed value of tx — t„iax ■ 

The solid curves indicate the evolution of t^v with fi^ 
which is given by equation (|19p. The two different curves 
correspond to the different assumptions about the value 
of the mode decay time due to non-linear effects, tni- The 
curve which is uppermost at small values of £ corresponds 
to the case a) where we assume that « t^v while the 
other solid curve corresponds to the case b) for which tni is 
much less than both tdec and t^v ■ One can see from Figures 
[2] and [3] that the behaviour of the characteristic time-scales 
is rather similar for these two cases. When £ is very small, 
rotation of the star is small, and both curves are close to 
the dot-dashed line. Also, for small values of £ we have 

tr = tev 

When £ increases, the star is spun up by tides and 
tr gets larger. When £ attains a value £gw ^ 10'* the 
solid curves, which for smaller values, correspond to tidally 
driven evolution cross the dotted line representing evolution 
controlled by the gravitational radiation time scale taw- 
For larger values of £, the orbital evolution is mainly de- 
termined by emission of gravitational waves, and we have 
tev = tew < tr- In this regime the star continues to spin up 
by tides, and the tidal time scales shown by dot dot dashed 
curves get larger with increasing of £. When the dot dot 
dashed curves cross the inclined dashed line, the binding 
energies are equal to Edis given by equations ([M)) and (pS)) . 
and, according to our criterion, the star is disrupted at en- 
ergies E ^ Edis- 

In Figure 2 we show the same quantities calculated for 
a star with a smaller value of orbital angular momentum 
corresponding to 77 = 4. In this case the evolution is always 
controlled by tidal effects and tev = tr for the whole range 
of energies shown £ < £dis and, accordingly, the the dot 
dot dashed curves coincide with the solid curves. 

As follows from this discussion and Figures [2] and [31 the 
white dwarf has the possibility of surviving the tidal evo- 
lution when the energy scale corresponding to the onset of 
the stochastic instability Est — GM / {2ast) (see equation 
(|13p ) is smaller than Edis - the energy corresponding to 
break-up rotation. When Est < Edis, disruption of the star 
is avoided during the phase of stochastic instability. In this 
case, when E > Est the orbital evolution proceeds mainly 
through emission of gravitational waves. Therefore, for our 
purposes it is very important to establish under what con- 
ditions the inequality Est < Edis holds. 

As can be seen from Figures [2] and [3l there are four dif- 
ferent possibilities: case al) where t„; « tev and tides de- 
termine the orbital evolution when the orbital energy of the 
star is close to Edis] case a2) where tni ~ tev and the gravi- 
tational radiation determines the orbital evolution; case bl) 



when tni ^ tev and tides determine the evolution, and the 
case b2) where tni «C tev and gravitational radiation con- 
trols the evolution. Note that it is straightforward to see 
that case bl) is always associated with disruption of the 
star. Indeed, in this case the factor tx/tev entering equation 
([26)) is of order of unity. That means that the energy Edis is 
always small when compared to Est , for typical parameters 
of the system. Therefore, the condition Est < Edis cannot 
be satisfied in this case and the case bl) is not considered 
further. On the other hand, in the opposite case b2) the fac- 
tor tx /tev can be large when estimated at energies of the 
order of Edis, see Figure 1, and the condition Est < Edis 
can be fulfilled, see Appendix for details. In Appendix it 
is also shown that circularisation is possible for the cases 
al) and a2) provided that the orbital angular momentum 
of the star and, accordingly, 77, is sufficiently large. Thus, 
circularisation can, in principal, be achieved for the cases 
al), a2) and b2). 

The condition Est < Edis can also be reformulated in 
terms of semi-major axes. Using the expression (j23p we 
can easily find the characteristic semi-major axis adis cor- 
responding to the binding energy per unit mass Edis 

ad^s = 7^^ 5-1013 f ^"j ) X^,-'Ahcm.m 

'^Odis \tT J \tnl + tdec J 

As we discussed above when a < adis and the condi- 
tion for the stochastic instability ([T3|) holds, emission of 
gravitational waves/non linear effects cannot significantly 
reduce the energy and angular momentum of the mode 
and the mode amplitude build-up persists. In such a sit- 
uation a white dwarf is likely to undergo intensive tidal 
heating which may possibly lead to disruption of the star. 
Therefore, the condition derived from (j30p together with 
corresponding conditions for onset of the stochastic insta- 
bility will be used in order to discriminate between orbital 
evolutionary tracks leading to disruption of the star and 
the tracks eventually leading to formation of close quasi- 
circular orbits around the black hole. 

When the non linear dissipation time scale tni is small 
(the cases bl and b2), adis does not depend on tni and we 
have 

ad^s = 7^ = ^- 10^' f \n-^M,cm. (31) 

'iodis \tT J 

When tni ~ tev (the cases al and a2) we can obtain an 
explicit expression for adis with the help of (|2ip : 

ad,s = - 4.4 • \Q^^\^/^^Ji-Hdl'\4>^r^'^cm.. (32) 

The condition ast > adis (or Est < Edis) can be written 
as a condition for the parameter rj to be larger than a cer- 
tain value. For our discussion below it is more convenient 
to use the quantity 0(77) defined in equation ([8|) instead of 
77. Depending on which case is considered (i.e. al, a2 or b2), 
the inequality ast > adis can be written as a condition that 
the quantity 0(77) is larger than its values ipi, (f>2 and 03 
corresponding to the equality ast = a-dis for the respective 
cases. These values are explicitly obtained in the Appendix, 
see equation (jA.Sp for the case al, equation (jA.7p for the 
case a2 and equation (jA.lOp for the case b2. 
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6. The rate of tidal capture by the black hole 

In order to estimate the rate of tidal capture by the black 
hole we assume that the stars in the cusp are slowly rotat- 
ing. Therefore, at the beginning of orbital circularisation 
the effects of rotation are not important so we neglect it 
when estimating the rate of capture into circularising or- 
bits, setting \[' = 1 in the expression (|2T|) for the time scale 
of orbital evolution. Then from the discussion in section 5 
and data plotted in Figures [5] and [3] it follows that during 
this phase orbital energy loss resulting from tides is more 
important than that due to gravitational radiation. 

The rate of capture into circularising orbits, as a func- 
tion of initial binding energy per unit mass of the star, 
depends on whether the tidal time scale ((2T|) is larger than 
the orbital period Porb{£)- Equating (pij) to ^ we find the 
binding energy per unit mass £i where these two times are 
equal to be given by 



£:i = o.7r'W( 



and we recall the relation between 77 and the O = J / Jc, 
that applies when the condition M < Merit as given by 
pl|l is fulfilled is e : 



1/3 



6.1. Balance of time scales for tidal evolution and orbital 
evolution due to distant encounters 

At first we assume that £ > £1 for a range of angular mo- 
menta of interest. A star gets tidally captured by the black 
hole when the characteristic time scale of evolution of its 
angular momentum due to distant gravitational encounters 
with other stars in the cluster 



J2 



orb-) 



(34) 



is larger than the tidal time scale tr (recall that the dis- 
persion of angular momentum j| is given by equation (O). 
Equating tj and tx we find the size of the " tidal circulari- 
sation" loss cone in angular momentum space as a function 
of the dimensionless binding energy £, 



Jt.( 



tr 



orb 



(35) 



Note that equation ([35]) defines Jt.c. implicitly since this 
quantity enters the right hand side through the dependence 
of tx on rj. To deal with this we define Qt.c. — Jt.c/ Jt, 
where we recall that Jt = \/2GMrT. Then we have 
©T.c. = V^^'^: see equations © and pU)l . Using the above, 
an equation for Qt.c. may be obtained from equations (O, 
(120]) and ^ in the form 



Qt.c. = 1.14Ai-i/3Mf ^Ai/20i/2(e3,c. 
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(36) 



where A* — A2/8.5. It follows from equation ((20|) that in 
order to find Qt.c. and hence the size of the tidal loss cone 
we need to evaluate the ratio AEt/E when J = Jt.c- 
This quantity grows from small values as the orbital period 
increases from small values, and accordingly the binding 
energy decreases, until it reaches unity when tj = Port- At 
that point equation (|33p applies as well as the condition 
that 



Jr 



T.C. 



(Aj2 



(37) 



These two equations define a critical energy Ecrit- 

When E < Ecrit the tidal evolution time scale is ap- 
proximately equal to the orbital period and therefore the 
size of the tidal circularisation loss cone is determined by 
equation (j33p which is, in this case, considered as an im- 
plicit equation for Q = Jt.c./ Jcone for a given value of E. 
Stars with these energies which have J < Jt.c. calculated 
in this way have relatively large energy changes induced 
within a single orbit. 

When E > Ecrit the number of stars tidally captured 
by the black hole per unit of time, with typical energies of 
order of some E, Nt.c.{E), can be estimated to be 



NT.C.iE) - KdiE, JT.cKEJT.c/i^ltT), 



(38) 



where we assume that the distribution function of the white 
dwarfs is given by equation Using equations (O 
and ([55]) we see that, in this case. 



(39) 



(33) Nt.c. « {Aj2fE cc E~^'^ 



is mainly determined by stars with small energies. Note 
that we take into account the fact that Ai depends only 
logarithmically on the binding energy when deriving the 
relation see equation ([55)1 . 

For the case E < Ecrit it can be easily seen that 
tj<tT~ Porb for any value of J. This situation is analo- 
gous to the full loss cone regime in the theory of tidally dis- 
rupted stars (e.g. Frank & Rees 1976, Lightman & Shapiro 
1977). In this case the rate of tidally captured stars can be 
estimated as 



Nt.c.{E) ^ N^aiE, JT.c)EJT.c/{^iPorb) « E, 



(40) 



where here we take into account the fact that the size of loss 
cone, now implicitly determined by equation ([33]), depends 
logarithmically and hence slowly on E. 

6.2. Determination of the critical binding energy per unit 
mass 

Thus, the rate of tidal circularisation is mainly determined 
by energies E Ecrit (Novikov, Pethick & Polnarev 1992), 
and therefore it important to obtain an explicit expression 
for the value of Ecrit- As we discussed above this value is 
determined by the condition that the timescale for orbital 
evolution due to tides be equal to the orbital period while 
being on the boundary of the tidal loss cone. It is deter- 
mined by equations (|33p and ([37]) which can be rewritten 
in the form 



£crit — 



GM 



: Ecrit —3-10 



where x = /iri/(AAf4) and 4'crit = 
equation ^ after setting 77 = Q^ fj , together with 



0(63 



T.C 



(41) 

) is found from 



^T.C. 



297r (ra 



£. 



-5/2 



crit 



(42) 



* Note an alternative form of equation (|38|l has already been 
provided in other studies, see e.g. Hopman & Alexander 2005, 
their equation (17) and references therein. To obtain this we 
write Nt.c.{E) ~ N^^iE, Jt.c)EJt.c / (trAi) = Ntot/itrAi), 
where Ntot = N^^{E, Jcirc)JcircE characterises the total num- 
ber of stars having binding energies per unit mass of or- 
der of E, we introduce the global relaxation time scale tr = 
J^i^^Porb/ {^j2Y and use equations dH) and (f551) . 
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Substituting Ecrit from (|¥T|) to we obtain a single equa- 
tion for rj and thence £crit in the form 



where 



(43) 



(44) 



and A* = A2/8.5. Using the exphcit expression for (j){r]) and 
the relation rj = Qj, q we rewrite (|43p in another form 



r/w 3.72 + 0.462 log 77 + 0.219 logy. 



(45) 



An approximate solution to equation (|45p valid for the 
range of interest of the parameter y ~ 0.05 — 10 with accu- 
racy of order of 10 per cent is 

ry« 4.3 + 0.22 logy, (46) 
and 

</.,„t(r;)«2-10V/'- (47) 

As seen from equation (|T7|) the size of the tidal circular- 
isation loss cone given by Qt.c. = V^^^ corresponding to 
the critical energy Ecrit is larger than the size of the tidal 
disruption loss cone given by O = 1. 

This is different from the situation where only the effect 
of emission of gravitational waves is taken into account as 
a mechanism for changing the orbital energy. In that case 
the energy scale analogous Ecrit may also be obtained from 
the requirement that the orbital evolution time t^v — taw 
be equal to the orbital period on the boundary of the corre- 
sponding circularisation loss cone. Obviously, our estimate 
of the circularisation rate is valid only when the size of the 
circularisation loss cone is larger than the size of 'true' loss 
cone, Jcone corresponding to direct capture or disruption of 
stars by the black hole. However, for typical parameters ex- 
pected for gravitational waves to be the dominant process 
of orbital evolution {M ~ IO^Mq and ri ^ 1) the energy 
scale Ecrit is smaller than another energy scale Econe de- 
fined by the condition that the sizes of the circularisation 
and true loss cones coincide at E = Econe- When E > Econe 
the stars are mainly destroyed by the black hole and the 
process of circularisation is strongly suppressed. In such a 
situation the capture rate is mainly determined by energies 
E ~ Econe (Hopman & Alexander 2005). Because of this 
difference, the tidal circularisation rate in our case, (see 
equation ([STjl below) differs from the estimate of the cap- 
ture rate due to emission of gravitational waves provided 
by Hopman & Alexander 2005, see also Novikov, Pethick 
& Polnarev 1992. 

Equation (|¥T|l gives 



150 



3/5 



We have explicitly 

£cru^l50Al/'{^y^%'^''M-'/'\ 



(48) 



(49) 



It is important to note that the critical energy Scrit ^ 1 
for our problem and therefore, the circularisation rate is 
mainly determined by scales well within the cusp around 
the black hole. Also, equation ((46|) tells that ~ 4 as we 
have mentioned above. 



6.3. Rate of capture into circularising orbits 

The maximal rate of capture , NJp^^ can be estimated a 
little more accurately as 



"T.C 



cJt.c. jY 
'-Jcrit I dJ - 



Jaon 



orb 



(50) 



where N^.d is given by equation (jj]) with E = Ecrit and 
Jcone is defined in section 2.4. Let us recall that in our 
problem Jcone — Jt- Integrating ((50|) we get 



N^C = -^5q-\'^)£cr^tna, 
v2Ai Ta 

where 

F = (e2,.^.(ine|.c.-i) + i)/4«o.2, 



(51) 



(52) 



and we assume that Qt.c = V'^^^ ~ 4^/'^ in the last equality. 
In astrophysical units the expression (|5ip has the form 

iV^"^- ^ 2.5-10-«5o.iA^/'A8-V"'/'Mf ^'"A^/'Vo^f /i%r-i.(53) 

where Sq,i — 5/0.1 and Ag = Ai/8. Note that we use ro.i = 
ra/{0.1pc) — lOri in the last equation assuming that a 
typical size of the radius of influence in globular clusters is 
of the order of O.lpc. 



6.4. Rate of capture of stars that will circularise without 
disruption 

As discussed in section [51 a white dwarf has a possibility 
of surviving the process of circularisation, and accordingly, 
settling on a quasi-circular orbit only when a typical value 
of the parameter 77 (or equivalently the pericentre distance) 
is sufficiently large. The tidally captured white dwarfs with 
formation rate iV^g? have a typical value of rj given by 
equation (|46[) and, accordingly, (pcritiil) given by equation 

(El- 
According to the discussion in section [5] and the 
Appendix, such white dwarfs have a considerable survival 
probability, for cases al) and a2) corresponding to large tni, 
when 

4>crit > (Pmax = niax(0i, (^2), (54) 

where 0i and 02 are determined by equations (|A.5|1 and 
(IA.7[) . respectively. For case b2) corresponding to small tni, 
for survival we require that 



> 



(55) 



instead of where 4>3 is given by equation (jA.lOp . 

When the above conditions are not satisfied, the typ- 
ical binding energy of white dwarfs having a considerable 
probability of survival, Es , has to be larger than the criti- 
cal value obtained from equation and so their capture 
rate is reduced by a factor [Eg/ EcritY^"^ when compared to 
that given by equation ([55)1 (see equation ([55)) ). Accordingly 
their capture rate is given by 



N. 



T.C. 



^Txf.i^s/ Ecrit) 
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(56) 



The binding energy per unit mass. Eg, is related to the 
quantity (l)max through the 'loss cone condition' p6[) . 
Taking into account the fact that the dependence of Qt.c. 
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on Es and (t>max is logarithmic, we approximately have 
Es^^ oc (f>max, and therefore equation (|56p can be written 
as 

Nt.c. « N^.ai^cru/<l>^ax), (57) 

and accordingly, the rate of formation is suppressed by the 
factor (pcrit I (j^max which is less than unity when conditions 
and ([55)) are not satisfied. As indicated below that is 
likely to be the situation for cases of interest. 

We go on to obtain explicit expressions for the rate of 
capture of white dwarfs that may survive the circularisation 
process in that case. This rate will be referred to as 'the 
circularisation rate'. 

First let us consider the case a) for which tni ~ tev 
As follows from equation (jA.8|) 02 > 0i when M > 2 ■ 
IO^Mq. Since the systems containing sufficiently massive 
black holes are likely to provide the bulk of the sources 
of gravitational radiation, we assume hereafter that M > 
2 • lO'^M© and set, accordingly, (j)max = 4'2- In this case the 
condition leads to 

r-o.i > SOA^ V"^^^! ■^ (58) 

where we use equations (jA.7|) and In the systems of 

interest this condition is not satisfied, and indeed the cir- 
cularisation rate is suppressed. 

We use equations (gTl), ^ and ^ to obtain the 
the suppressed capture or circularisation rate in the form 

N^^c. ~ 2.3 • 10-9<5o.iA^'-'AfO-57A-0-6Vo:}-5yr-i, (59) 

where hereafter we set the logarithmic factors A, = Ag = 1 . 
Substituting values of A = 1.4 and /i — 0.6 appropriate for 
a 'low density' white dwarf in (|59p . we obtain 

TVf c. - 8.2 • 10-i%iM4" "ro:i-^yr-\ (60) 

and in the case of the 'high density' white dwarf with A = 
0.84 and /i = 1, we obtain 

c. « 2.6 • lQ-^So.iM^'''ro l-^yr-\ (61) 

Now let us consider case b2) for which tni <C tev,tdec- 
In this case the condition ([55)1 can be expressed as 

ro.i > 10AV~'''/''Mi^', (62) 

where we have made use of equations (|A.10|) and (|i7)) . 
According to equation ((62)) the circularisation rate is equal 
to the rate of tidal capture ([55)1 when the cusp size is suf- 
ficiently large. For large cusp sizes the circularisation rate 
becomes small due to a strong decrease of the capture rate 
as it increases ( see equation (|55)) ). Thus such systems are 
not likely to contribute significantly to the total number of 
gravitational wave sources. Therefore, we assume hereafter 
that the inequality ((^ does not apply and multiply 
by the suppression factor {4>crit/4'3) in order to obtain the 
circularisation rate. Thus we have 

TVf c. ~ 6.2 • 10-9(5o.iAi^''/''Mf /3°A-6/5ro:i-^yr-i. (63) 

Substituting values appropriate to the two types of white 
dwarf in ((63)) . we obtain for the 'low density' case 

N^^c. ~ 1-8 • lO-''So,iMl%l-'yr-\ (64) 
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iVf c. ~ 7.7 • 10-^do.iMl-%.i^yr-' (65) 

for the 'high density case'. 

A comparison of equation ((55)) with equations ((5^ , ((50)) , 
((CT)) . ((55)) . ((M)) . ((55)) shows that our criterion of survival 
of the white dwarfs during the process of circularisation 
typically results in an order of magnitude decrease of the 
circularisation rate. However note that for case b2) when 
tni is small, the suppression is not significant for the 'high 
density' white dwarfs.. 

7. Probability of detection of a source of 
gravitational radiation 

The probability of formation of observable sources of grav- 
itational waves is determined by several important factors 
such as the rate of formation of white dwarfs on circular- 
ising orbits calculated above, the subsequent orbital evolu- 
tion of the white dwarfs, the abundance of globular clusters 
having sufficiently massive black holes and the properties of 
a gravitational wave antenna receiving the signal. In what 
follows we make a rough estimate of the probability assum- 
ing that the gravitational wave antenna has characteristics 
close to what has been indicated for the future LISA space 
borne gravitational wave antenna. 

LISA has its maximal sensitivity in the frequency range 
^ 10~^Hz for which the dimensionless amplitude of the 
gravitational waves, h, can be as small as 

h « 10-^^/^, (66) 

where Tobs is the time of observation, see e.g. Grishchuk, 
Lipunov, Postnov, Prokhorov and Sathyaprakash 2001 and 
references therein. Therefore, for simplicity, we assume that 
the white dwarf has an orbital period Porb ~ 100s when 
emitting the gravitational waves to be received by the an- 
tenna. The time Tobs can be smaller than or of the same 
order as the orbital evolution time T^y calculated for an 
orbital period 100s, and we assume hereafter that these 
are equal to each other, Tohs = T^v In order to estimate 
Tev we need to specify a typical value of eccentricity of the 
orbit, e, which, in turn, is determined by whether the or- 
bital evolution is governed by tides or by the emission of 
the gravitational waves. For such small orbital periods the 
evolution is likely to be determined by the latter effect and 
we can use the results obtained by Peters & Mathews 1963 
and Peters 1964 when estimating e. To do that we proceed 
as follows. First we estimate the orbital semi-major axis for 
an orbit with Porb = 100s from Kepler's law as 

a « 4:.7^^/^X-\t- (67) 

Next we use the dependence of the orbital angular momen- 
tum on e obtained by Peters, 1964 

J = OMH{e)e'^^^'^r)l/^JT, (68) 

where H{e) = (1 + i§e2)435/2299_ Equation §^ imphes 
that when e ~ 1 the orbital angular momentum is equal to 

1 /3 

J = r/f' Jt- Since a typical value of i] allowing survival dur- 
ing the process of tidal circularisation, according to section 
5 and Appendix, is slightly larger than the typical value 
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4, we set r/t — 5 hereafter. Now we can equate the 

to what is given 
5 obtain 



standard expression J — ^yGMa{^ — e 
in equation (|68p , and using equation (|67p with rjt 
an implicit equation for the eccentricity e 



,12/19 



(69) 



From this equation wc obtain e 0.56 and e 0.31 for the 
high/low density white dwarfs, respectively. The evolution 
time can be expressed in the form 



-2/3 p8/3 



(70) 



where P2 = Porfj/lOOs and the function F{e) specifies 
the dependence of T^v on e, see Peters 1964. We have 
F{e = 0) = 1, F{e = 0.31) w 0.8 and F{e = 0.56) « 0.2. 
Note, that an analogous expression for Te^, given in Ivanov 
2002 in his equation (79) contains an error in last expression 
on the right hand side and must be multiplied on additional 
factor M/(1O^M0). As a result, although his estimate of 
a probability of finding a source of gravitational waves in 
galactic centres containing black holes of masses ~ IO^Mq 
(equation (80) of his paper) is numerically correct for the 
parameters he uses, it has a different scaling with the black 
hole mass. It must be multiplied by an additional factor 
M/(1O^M0) which results in smaller values of the proba- 
bility for black hole masses in the range ~ 10** — IG^M©. 

Now, substituting (j70p in (|66p we obtain minimal values 
of the amplitude h, hmm, which can be detected by LISA 



/i^'^. «io- 

"■mm ~ '-^ 



where h^jl^^^ stands for the minimal amplitudes for the cases 
of low/high density white dwarfs, respectively, and we set 
P2 = 1. 

In order to find the total volume of space, available for 
detection of the signals by LISA, we should relate the quan- 
tities h^^^^ to the distance from a source, Robs- For that the 
amplitude h of radiation emitted by the source should be 
known as a function of distance and of parameters of the 
source. The appropriate relation has been obtained, e.g., 
by Nelemans, Yungelson and Portegies Zwart 2001. In our 
units it has the form 



5 • 10'^*M, 



24 ,.1/3 



(71) 



h{Robs) « 5 • 10-''fif{n,e)M^'-'P 



2/3 p-2/3 IMpc 



(72) 



where the function f{n, e) contains information about the 
amplitudes associated with harmonics of different order n. 
It is related to the function g{n, e) calculated by Peters & 
Mathews 1963 through /(n, e) = ^ g{n, e)/n 0. We have 
f{n, e — 0.56) « 0.2 for several harmonics with n = 1,2,3 
and f{n = 2,e — 0.31) w 0.4. Now we can equate h{Robs) 



to h. 



(1,2) 



and find the maximal distance to a source which 



can be, in principal, detected by LISA 



R 



(0 

obs 



1200M|/^Mpc, R);ii = 200M;'"Mpc, 



(h) 



ri/3, 



where the upper indices {I, h) stand for the cases of 
low/high density white dwarfs hereafter. 



^ Note a misprint in Appendix of Peters & Mathews 1963 in 
their expression for g{n,e). The fifth term in the braces must 
contain the factor (4/e)^ contrary to the factor (4/e^)^ as in the 
text. 



In order to find the total number of potential sources, 
A^tot: we use an estimate obtained by Portegies Zwart & 
McMillan 2000 for the average number of globular clus- 
ters in the Universe, na.c. ~ iMpc^^ for the value of the 
Hubble constant H = 70km- -Mpc^^^ We obtain from 
equation ([73]) 



fna.c.AiRillf 



2 • 10^"AM4 



N^!:i « 108 AM4, 



(74) 



where A is the fraction of globular clusters having black 
holes with masses ~ M4. The total probabilities of finding 
of a source, 

P^(l:2) 

can be estimated as products of the 
circularisation rates given by equations (j60p and (|6ip for 
case a) and equations ([M]) and Ip5|) for case b), the total 
numbers of the sources available and the operational time 
of the LISA mission taken to be equal to 2yr. 

For case a) with tni ~ tev, using equations (pO)) . (pT|) 
and (|74p we have 



Pr('' w 30a, Pr'-^'^ w 0.4a, 
where 

a = A5o,iM} 



'0.1 



(75) 



(76) 



It is seen from equation (j75p that the high density white 
dwarfs of Sirius B type give rather small probability of de- 
tection. On the other hand, the 'ordinary' white dwarfs 
could provide, in principal, a source of gravitational waves 
provided that our criterion of white dwarfs survival is valid 
and that 



a > 1/30. 



(77) 



It is interesting reformulate the criterion (|77p in term of a 
typical velocity of stars near the radius of influence of the 
black hole, a ~ ^ GMjva- We have from equations ((76)) 
and (173 



(78) 



^ 1 /'20km/s 



where we neglect a very slow dependence on the black hole 
mass, a{a,M) (x M°-°^. Of course, the black hole mass 
must be sufficiently large for the validity of the assumptions 
leading to equation ([77]) . 

For case b) with tni <^ tev , tdec we proceed in a similar 
manner but use equations (|64p and ([55)1 instead of equations 
(ISO)) and ((CT)) . respectively. We find 



Pr(') « 70a, Pr'-''^ 
where 



1.5a, 



a = A(5o.iM, 



(73) and 



A > 



2.2-1.5 
'0.1 • 



20fcm/s 



70So: 



(79) 



(80) 



(81) 



Note that the average number of globular clusters is mean- 
ingful only when -Robs > IQMpc, where 70Mpc is a typical dis- 
tance between the elements of the Large Scale Structure in the 
Universe. 
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instead of equation (|78p for a 'typical low density' white 
dwarf. 

As for case a) 'typical' white dwarfs have a much larger 
probability of detection. However, for case b) this proba- 
bihty depends explicitly on the black hole mass when the 
cusp size is expresses in terms of the typical stellar velocity 
dispersion a. As follows from equations ([75]) , (|80p and (|5T|) 
when M4 ~ 1 and a ~ 2Qkm/s this probability is about 
two times larger than that corresponding to case a). But 
the probability ratio decreases with the black hole mass if 
cr is kept fixed. 

8. Discussion 

In this paper we have given a qualitative analysis of the 
problem of tidal circularisation of white dwarfs in globu- 
lar clusters containing black holes with masses in the range 
10'^ — 10'* M0 and estimated the rate of production of white 
dwarfs that begin to circularise their orbits for specified 
parameters for the cluster, white dwarf and black hole, see 
equation (|53p . We also proposed a simple criterion for 'sur- 
vival' of a white dwarf during the stage of orbital circulari- 
sation, and found production rates for white dwarfs, which 
can, in principal, settle down in a quasi-circular orbit thus 
forming sources of gravitational radiation, see equations 
(|59II61I) and (|63ll65p . These rates are order of magnitude 
smaller than that given by equation (|53p . We made a simple 
estimate of the probability of detection of systems contain- 
ing white dwarfs on quasi-circular orbits and found that 
these systems can, in principle, be detected by LISA pro- 
vided that the globular clusters containing black holes are 
sufficiently abundant and their stellar velocity dispersion 
near the radius of influence of the black hole is sufficiently 
large, see equation ([78|) and ([8T|) . 

Our results should be treated with caution. Many pro- 
cesses occurring during the formation of circularising stars 
and during ongoing tidal circularisation depend strongly on 
the parameters of the problem and this leads to uncertain- 
ties in the results. Let us consider some of them. 

1) In a quantitative approach a realistic model of a white 
dwarf with non-zero temperature should be used. The 
eigen spectrum of perturbations should include the grav- 
ity/inertial modes. The realistic model should also take 
into account different channels of mode decay and a self- 
consistent change of structure of the white dwarf as a result 
of radiative cooling/heating due to dissipation of the mode 
energy. 

2) Our results indicate that the stars can achieve very high 
rotation rates during the process of ongoing tidal circular- 
isation. The possibility of survival of white dwarfs during 
this stage critically depends on tidal exchange of energy 
between the orbit and the modes of a fast rotating star. It 
will be very important to generalise well known results of 
the theory of dynamic tides in non-rotating/slowly rotating 
stars to the case of fast rotators. 

3) As follows from this discussion, our results depend on 
the rate of mode energy and angular momentum transfer 
to the star. Possible channels for this transfer, especially 
through dissipation resulting from non linear effects should 
be further investigated. 

4) The dependence of the circularisation rate on the cusp 
size and black hole mass should be checked with numerical 
N-body simulations. Also, a more realistic distribution of 
the white dwarfs over the cluster as well as a change of 



their relative fraction due to stellar evolution should be 
accounted for in a numerical study. 

5) We have adopted a rather simplistic criterion for de- 
tection of the gravitational waves produced by the white 
dwarf orbiting close to a black hole, assuming that the am- 
plitude of the signal is equal to the noise amplitude of the 
detector. Higher S/N detection threshold would result in a 
much smaller probability of detection. On the other hand, 
a more sophisticated strategy of signal processing could im- 
prove our estimates. 

We have shown that during the first stage of orbital 
circularisation the orbital energy transferred to the funda- 
mental mode of stellar pulsations may be radiated away by 
gravitational waves. That means that the globular clusters 
may contain sources of gravitational radiation with a typ- 
ical frequency of the order of the eigen frequency of the 
fundamental quadrupole mode, of the order of a few Hz. 
The possibility of detection of such sources needs further 
investigation. 

Finally we comment that in a recent paper Baumgardt, 
Hopman, Portegies Zwart, and Makino 2006 have calcu- 
lated the tidal capture rates for stars of different types with 
help of N-body simulations, during first 12Myrs of evolu- 
tion of the system. They have considered black holes with 
masses in the range ~ 10'^ — 4 • IO^Mq. Giant stars and low 
mass main sequence stars with mass < 0.4Mq have been 
modelled as n = 1.5 polytropes. The low mass stars seem 
to be absent in the cusp due to the effect of mass segre- 
gation operating during formation of the cusp. However, 
their capture rate of red giants can be compared with what 
is given by equation (|53p . They have obtained a capture 
rate for red giants of the order of a few events per run. 
Assuming that a typical red giant has a radius '-^ lO^i?©, 
mass ~ IOA/0, taking the black hole mass and the cusp 
size to be equal to 3 • IQ^Mq and to O.lpc, respectively, 
and substituting these values into equation (j53p , we obtain 

Nr.g. ^ 3 • 10~^ S R, G. yr~^ , where Sr g. is the number frac- 
tion of the red giants in the cusp. Our results appears to 
be in full agreement with the results of N-body simulations 
provided that (J^.g. ^ 1/100. 

Acknowledgements 

We are grateful to A. G. Polnarev for fruitful discussions 
to A. G. Doroshkevich for useful remarks and to the referee 
for valuable comments. PBI has been supported in part by 
RFBR grants 04-02-17444 and 07-02-00886. 

Appendix A: Conditions for the transition from 
stochastic instability 

As indicated in section [3771 there is a phase of evolution at 
large semi-major axis subsequent to tidal capture where the 
internal pulsation mode energy may increase in a stochastic 
manner. This requires the semi-major axis to exceed a^t 
above which the pulsation mode does not maintain phase 
coherence between successive pericentre passages. 

Here we obtain an explicit expression for the semi-major 
axis flsi defined through equation (|13p . At first let us assume 
that AEt > AEgw during the orbital evolution. This sit- 
uation is illustrated in Figure 2. In this case we substitute 
equation ^ in equation (I13p and obtain 

ast = ai]'> w 1.15 • 10"AAi"^/^Mf/^(</.*)2/5cm. (A.l) 
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Now let us consider the orbital evolution when A£'t < 
AEgw ■ This situation is illustrated in Figure 1 . In this case 
the semi-major corresponding to the onset of the stochastic 

(2) 

instability, Ust = dst can be obtained by substitution of 
equation (pS)) in equation (fT5|) with the result 

ai?' « 1.7 • 101^7714/15 AV^^'/^'Mycm. (A.2) 

Just after tidal capture occurs the orbital evolution is 
primarily due to tides. The orbit decays until a^t is reached 
but at this point it is possible that the orbital evolution is 
governed by gravitational radiation. As indicated in section 
[51 when the orbital semi-major axis is of the order of ast 
there are three non trivial possibilities: case al) where i„; ~ 
tev and tides determine the orbital evolution; case a2) where 
tni ~ tev and gravitational radiation determines the orbital 
evolution; and case b2) where t„i <C tev and gravitational 
radiation controls the orbital evolution. Each of these is 
considered in more detail below. 

A.l. The case al: tides dominate and tni is large 

As indicated in sections 4 and 5, when tni ^ t^v a white 
dwarf has a possibility of surviving the process of tidal cir- 
cularisation only if it is cooled efficiently by emission of 
gravitational waves for values of its semi-major axes > a^j^ . 
Therefore, in order to reach a quasi-circular orbit, the white 
dwarf must have orbital parameters such that a]^^ > adis 
where adis is given by equation (|32p . This inequality leads to 
an inequality of the form > 771 ( or equivalently the fixed 
pericentre distance must exceed a certain value), where 771 

is defined by the condition a^j^ — adis- In general, 771 is a 
monotonic function of the stellar rotation rate with smaller 
values of 771 corresponding to larger values of Sl^. This is 
because tides weaken with increasing £7^ enabling safe cir- 
cularisation starting with smaller pericentre distances. 

Thus in order to take into account all possible evolu- 
tionary tracks leading to formation of circular orbits, we 
calculate the lower boundary of 771 assuming that fi^ — 
Slfcr ~ 0.5il*. In this case the factor (^^E") entering (jA.ip 
takes the form 

{(j)-^) w - exp (3.6877 - 2.74). (A.3) 
'7 

Substituting equation (|A.3p in equation (jA.lD and equating 
((5^ and (|XT|) . we obtain 

7^1 « 4 + 0.27(ln?7i +lnZi), (A.4) 

where Zi = X^'' I ii-''^ I Ml^ ' \ 

In the subsequent discussion we use the quantity 0i = 
(^(771) defined in equation ([5]) rather than 771. This quantity 
can be approximately represented in the form 

« 2.4 • lO^Z?-^ « 2.4 • 103a1V"° ''M4-° i^5^ (A.5) 

with the difference between this power-law fit and the so- 
lution of equation (|A.4[) being less than or of the order of 
5 per cent for all interesting values of Z\ . 

A.2. The case a2: gravitational waves dominate and tni is 
large 

Now let us assume that gravitational waves determine the 
orbital evolution when a rapidly rotating white dwarf has 



its semi- major axis near to adis and tni ~ t^v ~ ^gWi see 
Figure 1. Proceeding as in the preceding section, we obtain 
an equation for the critical value of 772 that is analogous to 
equation (|A.4p in the form 

772 w 4.88 - 0.11 In 772 + 0.27 In Z2, (A.6) 

where Z2 = A/i"^/^ Af^^^. The corresponding value of 02 = 
0(772) can also be, with a high accuracy, represented in the 
power-law form 

02 « 5.4 • lO^zf « 5.4 • 103A°-SV"° ''M4"-26. (A.7) 

The value of the parameter 77 and thus 0(77) which di- 
vides, according to our criterion, orbits leading to a safe 
circularisation from orbits leading to a possible disrup- 
tion of the white dwarf, is obtained from the condition 
0(77) > 777,ax(0i, 02). We note that the condition for grav- 
itational radiation to dominate, 02 > 0i, leads to the ap- 
proximate condition on the black hole mass that 

A'U > A/, = 0.155 A, (A.8) 

where we have neglected a very weak dependence of 
on /i of the form Af, cx the latter quantity being 

replaced by unity. 

The inequality (|A.8p implies that for a sufficiently large 
values of the black hole mass and the rotation of the white 
dwarf the orbital evolution is mainly determined by emis- 
sion of gravitational waves when the semi-major of the orbit 
is close to the characteristic semi-major axis corresponding 
to the onset of stochastic instability. For the parameters 
corresponding to a 'typical' white dwarf with /i — 0.6 and 
A = 1.4, we have A/, « 0.22, and for a dense white dwarf of 
Sirius B type with /i = 1 and A — 0.84, we have Af* = 0.13. 
Both characteristic values of A/* are smaller than the black 
hole masses of systems expected to produce a large amount 
of gravitational radiation during the last stages of orbital 
evolution of the star. Therefore, the condition 77 > 772 is 
more important for our purposes. 

A.3. The case b2: gravitational waves dominate and tni is 
small 

Now let us consider the case when tni ^ t^v^tdec- The 
discussion proceeds in a similar way to the previous case 
but now adis is given by equation pip . As we have men- 
tioned above it can then be shown that gravitational waves 
dominate the orbital evolution when the semi-major axis is 
~ adis provided that the black hole mass exceeds ~ 10^ Mq, 
and accordingly, tev ~ taw- 

Therefore in order to obtain an equation for the param- 
eter rj = rj^ separating orbits leading to safe circularisation 
from those leading to disruption of the star, we equate ex- 
pressions ^ and (|A.2p . and use equations 0, dH]) 
and (|A.3p . to obtain 

7/3 = 3.5-|-0.651n7;3 -t-0.271nZ3, (A.9) 

where Z3 — X^/'^ M^^^^^ . An approximate solution to 
this equation can be found from 

03 = 0(773) « 2 • 103A«/5^-24/25m732/'^ (A.IO) 
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